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Abstract. We study the determinant det(/ — -Kpu) of an integrable Fredholm opera- 
tor Kpu acting on the interval (— s, s) whose kernel is constructed out of the ^-function 
associated with the Hastings-McLeod solution of the second Painleve equation. This 
Fredholm determinant describes the critical behavior of the eigenvalue gap probabilities 
of a random Hermitian matrix chosen from the Unitary Ensemble in the bulk double 
scaling limit near a quadratic zero of the limiting mean eigenvalue density. Using the 
Riemann-Hilbert method, we evaluate the large s-asymptotics of det(J — Kpu). 



1. Introduction and statement of the results 
Let M{n) be a unitary ensembles of random n x n Hermitian matrices M = {Mij\ 



<*S f M equipped with the probability measure, 

P^ N \M)dM = ce- NtrV ^dM, c f e- NtlV ^dM = l. (1.1) 

M(n) 

Here dM denotes the Haar measure on M.(n) — W 1 , N is a fixed integer and the 
potential V : R — > R is assumed to be real analytic satisfying the growth condition 

V(x) , \ 

- r — y oo as \x\ —> oo. (1-2) 

m(x z + 1) 

The principal object of the analysis of the model is the statistics of eigenvalues of the 
matrices from M (n) . A classical fact [8] , [25] is that the eigenvalues form a determinantal 
random point process with the kernel 

n-l 

K n , N (x,y) = e~fy(-) e -f v ^J2Pi( x )Pi(y)^ (1-3) 

i=0 

where Pj(x) are polynomials orthonormal with respect to the weight e ~ NV ( x ) f 



V 



(x) Pj (x)e- NVi - x) dx = dij, Pj (x) = K jX j + ... (1.4) 



In particular, one of the basic statistical characteristics, the gap probability, 

Eu,n(s) = Prob(M 6 Ai(n) has no eigenvalues in the interval (— s, s), s > 0) 
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is given by the formula, 

n . S S 

" ^ (—~\x> r r 

E n ,N(s) = ^2 — T~ / " / det { K n,N(xk,xi) j k l=1 )dx 1 ■ ■ ■ dxj 
;=o ^' 

— s —s 

= det(I-K ntN ), 

where K n ^ is the trace class operator acting on L 2 ((— s, s), dx) with kernel K h ,n(x, y). 

Assumptions (jl.2p on the potential V{x) ensure [12] (see also [8] for more on the 
history of the subject) that the mean eigenvalue density —K n ^{x,x) has a limit, 

lim — K n n(x, x) = pv(x) > 0, (1-5) 

n,N— >oo 71 ' 



whose support, = {x £ M. : py(x) > 0} , is a finite union of intervals (simultaneously, 
Pv( x ) defines the density of the equilibrium measure for the logarithmic potentials in 
the presence of the external potential V). The limiting density pv( x ) is determined by 
the potential V(x). At the same time, the local statistics of eigenvalues in the large n, N 
limit satisfies the so-called universality property, i.e. it is determined only by the local 
characteristics of the eigenvalue density py ([26], [S],[13]). For instance, let us choose a 
regular point x* S Sy, i-e. pv(x*) > 0. Then the bulk universality states that 

lim -^--K n Jx* + -^r-ryx* + -^—) = K sm (\, - ^ X ~ f (1.6) 
n->oo npvyx*) \ npv\x*) npy[xr) ) tt{\ — p) 



uniformly on compact subsets of M, which in turn implies [13] that for a regular point 

x , 

lim Probl M G M{n) has no eigenvalues G [x* 



V_ " " v ' & "V npy(x*Y npv(x*) 

= det(/-^ sin ), (1.7) 

where -K" s i n is the trace class operator on L 2 ((— s, s); dX) with kernel K s i n (\, p) given 
in (jl.6p . (This result was first obtained for GUE in the classical works of Gaudin and 
Dyson.) The Fredholm determinant in the right hand side of (jl.7p admits the following 
asymptotic representation [15], 



lndet(/-^ sin ) = -^-iln(vr S ) + ^ln2 + 3C'(-l)+0(^ 1 ), s -)• oo, (1.8) 

where C'(^) is the derivative of the Riemann zeta- function (a rigorous proof without 
the constant term was obtained independently by Widom and Suleimanov - see [TT] 
for more historical details; a rigorous proof including the constant terms was obtained 
independently in [21], |16] - see also [H]). This remarkable formula yields one of the most 
important results in random matrix theory, i.e. an explicit evaluation of the large gap 
probability. 

Equation (jl.7p shows that in double scaling limits the basic statistical properties of 
hermitian random matrices are still expressible in terms of Fredholm determinants. This 
is also true for the first critical case, when pv{x) vanishes quadratically at an interior 
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point x* G Sy. However, in this situation the scaling limit is more complicated [3], |6j. 
Let (x*) = p' v (x*) = 0, p'y(x*) > and n, N — > oo such that 



lim n 2 / 3 f--0 = C 

n,Af->oo \N ) 



exists with C £ R. Then the critical bulk universality guarantees existence of positive 
constants c and c\ such that 

U m — TTT^ra,^ ( x * H TH^^H TH ) = ^Pii(A,/i;x) (1.9) 

n,iv ->oo cra i / d \ era 1 / 13 en 1 ' 6 J 

uniformly on compact subsets of K where the variable x is the scaling parameter defined 
by the relation 

lim re 2 / 3 ( — — 1 ] = xc\. 
n,N-Hx \N J 

Here the limiting kernel Kpu(X,/i;x) is constructed out of the ^-function associated 
with a special solution of the second Painleve equation. The precise description of the 
kernel Kpn(\, fi; x) is as follows. 

Let u{x) be the Hastings-McLeod solution of the Painleve II equation |2Uj . i.e. the 
unique real- valued solution to the boundary value problem 



xu + 2u 3 , u(x) 



Ai(x), x — > +oo; 



x — > — oo, 



where Ai(x) is the Airy-function (the solution u(x) is in fact uniquely determined by 
its Airy-asymptotics at x = +oo). Viewing x, u = u(x) and u x = du(x)/dx as real 
parameters, consider the 2x2 system of linear ordinary differential equations, 

s. T . ,,,, N ,., 2 A -\ ( u\ f-ix-2iu 2 —2u x 



dx -A(x, x) *,A(x, x) = -iiy* 3+ «x{: u ^ _ 2 - jx ; 2 7„ 2 j.(i.io) 

Let ^(A) = ^(A, x) be the fundamental solution of system (|1.1U|) which is uniquely fixed 
by the asymptotic condition, 

m(\,x) = (l + 0(\- 1 ))e- i ^ xS+xX)a \ A^oo, 0<argA<7r. 



Then, the kernel Kpu(\, /i; x) is given by the formula, 

JTpnCA, m x) ee K PII (X, n) = —( MX,x)Mp,x)- ^,x)MX,x) \ 

2ir \ A — ji J 

where ?/>ii(A, x) and ^21 (A, x) are the entries of the matrix valued function ty(\,x) = 

{lpjk(X,x)}j,k=l,2- 

Remark 1. The function \I/(A,x) can be alternatively defined as a solution of a certain 
matrix oscillatory Riemann-Hilbert problem. The exact formulation of this Riemann- 
Hilbert problem is given in the next section. 

In this paper we study the Fredholm determinant 

det(/-K P „), (1.12) 
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where -fTpn is the trace class operator on L 2 ((— s, s); dX) with kernel (jl.lip . In virtue of 
(jl.9j) . this determinant replaces the sine - kernel determinant in the description of the 
gap-probability near the critical point x* , i.e. instead of (|1.7p one has that 



limProbl M G M.(n) has no eigenvalues E [x 



:,x + 



\ v cn 1 / 3 ' cn 1 / 3 . 

= det(/-tf Pn ), (1.13) 

as n, N — > oo and 

lim n 2 / 3 [ — — 1 ) = xc\. 

n,N^oo \N J 

(The proof can obtained in a same manner as the proof of the similar equation (21) in [10] 
with the help of the proper estimates from [4].) The main result of the present paper is 
the following analogue of the Dyson formula (jl.8p for the Painleve II determinant (jl,12p . 

Theorem 1. Let Kpu denote the trace class operator on L 2 ((— s, s); dX) with kernel 
(jl.lip . Then as s —> oo the Fredholm determinant det(I — -ftTpn) behaves as 

oo 

2 13/" 
lndet(7 - ifpn) = ~ s ^ x ~ tj^) 2 ~ ^ ms + (y ~ x)u 2 {y)dy 

X 

_I ln2 + 3C / (-l) + 0( S - 1 ), (1.14) 

and the error term in (|1.14p is uniform on any compact subset of the set 

{i£R: -oo<x<oo}. (1-15) 

We bring the reader's attention to the following two interesting aspects of formula 
(|1.14p . One is the appearance in the asymptotics of the Tracy- Widom function, 

F TW (x) = e -^^- x)u2{v)dy . 

The second important feature of the estimate (|l,14p is related to the For r ester- Chen- 
Eriksen- Tracy conjecture ([E], [5]; see also [2]) concerning the behavior of the large 
gap probabilities. The conjecture states that the probability E(s) of emptiness of the 
(properly scaled) interval (x* — s, x* + s) around the point x* satisfies the estimate, 

E{s) ~exp(-Cs 2K+2 Y (1.16) 



if the mean density p(x) behaves as p ~ (x — x*) K . This conjecture is supported by the 
classical results concerning the regular bulk point (k = 0, the sine-kernel determinant) 
and regular edge point (k = 1/2, the Airy- kernel determinant). For higher order critical 
edge points (k = 21 + 1/2, the higher Painleve I - kernel determinants), estimate (|1.16D 
follows from the asymptotic results of [Tj- Our asymptotic equation (jl . 14|) supports 
the Forrester-Chen-Eriksen- Tracy conjecture for the first critical case in the bulk, when 
k = 2. 

The proof of Theorem [T] is based on a Riemann-Hilbert approach. This approach 
(compare [21], |11| ) uses the given integrable form of the Fredholm operator, allowing 
us to connect the resolvent kernel to the solution of a Riemann-Hilbert problem. The 
latter can be analysed rigorously via the Deift-Zhou nonlinear steepest descent method. 
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Let us finish this introduction with a brief outline for the rest of the paper. Section [2] 
starts with the more details concerning the definition of the kernel (jl.lip followed by a 
short review of the Riemann-Hilbert approach for the asymptotics of integrable Predholm 
operators. We then apply the general framework to the Fredholm determinant det(J — 
Kpn) and formulate the associated "master" Riemann-Hilbert problem (RHP). We will 
also evaluate logarithmic s and x derivatives of the determinant det(J — K s ) in terms 
of the solution of this RHP and outline a derivation of an integrable system whose tau- 
function is represented by det(/— -Kpn) . Since the main objective of this paper is the large 
s asymptotics of det(J— -Kpn), we shall postpone a thorough discussion of the differential 
equations describing this determinant until the next publication. In sections 141 71 — TT2l 
following the Deift-Zhou scheme, we construct the asymptotic solution of the master 
Riemann-Hilbert problem. Comparing to the more usual cases, an extra "undressing" 
step is needed to overcome the transcendentality of the kernel Kpu(X, fi; x). Here, a 
crucial role is played by the aforementioned alternative Riemann-Hilbert definition of 
the function ^(A, x). The situation is similar to the one dealt with in [7]. The calculations 
of sections [13] and provide us with the asymptotics of lndet(I — Kpn) announced in 
(|1.14p but up to the constant term. In order to determine the latter, we will, in section 
[T5l go back to equation (jl.lip and look at the behavior of the kernel Kpu(X, /j,;x) as 
x —7- oo. We will see that in the large x limit, the kernel Kpu(X, /i) is replaced by the 
following cubic generalization of the sine kernel, 

, s - , N sin (|(A 3 - ^ 3 ) + x(A - ^)) , 
Kpn(X^) m- A- csin (A,^) = — ^ n / — > (1-17) 

We will then introduce a parameter t £ [0, 1] 

sin(|t(A 3 -/i 3 ) + x(A-/i)) 



K csin (X,fi) i-)- K csin (X,n) 



tt(A - n) 



and compute the large s behavior of det(7 — K cs [ n ) using again the Riemann-Hilbert 
approach. This will be done in sections [T71I24I This analysis will indeed produce the 
constant term in (|1.14p since det(I — K csm ) | reduces to the sine kernel with known 
asymptotics see (|1.8p . 



lndet(/ - K csin )\ t=0 = - j ln(sx) + ^ In 2 + 3C'(-1) + O^" 1 ), s oo 

uniformly on any compact subset of (|1.15|) . 

Remark 2. We do not address in this paper the question of the higher corrections to 
After the leading and constant terms are determined, the higher corrections can 
be in principal obtained by iterating the final R- Riemann-Hilbert problem (see section 
\14\ )- Alternatively, one can use the related to the determinant det(/ — Kpu) differential 
system which we have mentioned above, and which we intend to discuss in detail in our 
next publication. 

The asymptotic of the determinant det(J — i^csin) is of interest on its own. We show, 
that it is given by a formula ignoring the Tracy- Widom term in (|1.14p . That is, our 
second asymptotic result is the following. 
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Theorem 2. Let K cs - m denote the trace class operator on L 2 ((— s, s); dXj with kernel 
(jl,17|) . Then as s — )■ 00 the Fredholm determinant det(J — i^csm) behaves as 

lndet(J - K csin ) = ~s 6 - s 4 x - ^{sxf - jjlns - i In 2 + 3C'(-1) + O^" 1 ). (1.18) 
and i/ie error term in (|1.18p is uniform on any compact subset of the set f)1.15[) . 

2. RlEMANN-HlLBERT APPROACH - SETUP AND REVIEW 

The classical theory of ordinary differential equations in the complex plane implies 
that system (|1.1U|) has precisely one irregular singular point of Poincare rank 3 at infinity. 
This observation leads to the existence of seven canonical solutions ^(A) which are fixed 
uniquely by their asymptotics (for more detail see [17j ) 

tf n (A) ~ (l + 0(A~ 1 )) e - i (f A3 +- A )-3 5 a _> A e o n 

where the canonical sectors O n (compare Figure [1]) are defined by 

O n = {A e C I arg A G (|(n - 2), |n) ,n = 1, . . . , 7}. 




Figure 1. Canonical sectors of system (|1.10p with the dashed lines in- 
dicating where Re A 3 = 

Moreover the presence of an irregular singularity gives us a non trivial Stokes phe- 
nomenon described by the Stokes matrices S n : 

S n = 0„(A))-V„ +1 (A) 
In the given situation (jl.lip (see again [17] ) these multipliers are 

Si = jj , S2 = jj , S4 = ^0 l) ' ^ 3 = ^ ^ 5 = ^ 6 = 
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hence the required solution in (jl.lOp is the second and third canonical solution ^(A, x) = 
^2(A,x) = \I/3(A, x) with asymptotics 

*(A,x) ~ (l + 0(X~ 1 ))e- i{ l x3+xX)a3 , A^oo, < arg A < vr (2.2) 



and Stokes matrices as in (|2.ip . Now that we have defined the integral kernel let 
us locate its structure within the algebra of integrable Fredholm operators and discuss 
the underlying Riemann-Hilbert approach. The given kernel belongs to an algebra of 
integrable operators first introduced in [21 j : Let T be an oriented contour in the complex 
plane C such as a Jordan curve. We are interested in operators of the form XI + K on 
L 2 (r), where K denotes an integral operator with kernel 

K(X^) = E ^^ i{X)ki ^\ V/,(AMA) = 0, MeZ>! (2.3) 
A " P i=i 

with functions fi,hi which are smooth up to the boundary of T. Given two operators 
XI + K , XI + K of this type, the composition (XI + K)(Xl + K) is again of the same 
form, hence we have a ring. Moreover let K t denote the real adjoint of K, i.e. 

A — f^i 

Our results are based on the following facts of the theory of integrable operators (see e.g. 
|11|). First an algebraic Lemma, showing that the resolvent of I — K is again integrable. 



Lemma 1. Given an operator I — K on L 2 (T) in the previous algebra with kernel 
Suppose the inverse (I — K)" 1 exists, then I + R = (I — K)~ l lies again in the same 
algebra with 

R{x ) = S£ifl(A)iW Y, m)Hi{X) = o (2 .4) 

and the functions Fi,Hi are given by 

Fi(X) = ((/ - i0~7*) (A), Ht(X) = ((/ - (A). (2.5) 

Secondly an analytical Lemma, which connects integrable operators to a Riemann- 
Hilbert problem. 

Lemma 2. Let K be of integrable type such that (I — K)~ l exists and let Y = Y(z) 
denote the unique solution of the following N x N Riemann-Hilbert problem (RHP) 

• Y(z) is analytic for z G C\r 

• On the contour T, the boundary values of the function Y(z) satisfy the jump 
relation 

Y + (z) = Y_(z)(l-2nif(z)h t (z)), z€T 

where f(z) = (fi(z), . . . , /at(z))* and similarly h(z) = (hi(z), . . . , h]\r(z)Y 

• At an endpoint of the contour T, Y(z) has no more than a logarithmic singularity 

• As z — > oo 

Y(z) = I + 0(z- 1 ) 
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Then Y{z) determines the resolvent kernel via 

F(z) = Y(z)f(z), H(z) = (Y\z)Y 1 h{z) (2.6) 

and conversely the solution of the above RHP is expressible in terms of the function F(z) 
using the Cauchy integral 

Y{z) = I - [ F(w)h\w)-^-. (2.7) 
J w - z 

r 

Let us use this general setup in the given situation (jl.lip . We have 

KpuM = tmE>, ' HM-'(*f<»'\) (2.8) 

A-/i V2vr V^2i(A,x)y V2ir\-tpu{fi,x)J 

hence the Y-RHP of Lemma 2 reads as 

• Y(A) is analytic for A £ C\[— s, s] 

• Orienting the line segment [— s, s] from left to right, the following jump holds 

Y + (\) = Y-(\)( 1 + %Wf^ V \e[-s,s] 

• At the endpoints A = ±s, Y(A) has logarithmic singularities, i.e. 

Y(X) = 0(ln(A=F*)), A -> ±a 

• As A —7- oo we have 

y(A) = /+^i + 0(A" 2 ). 

A 

The given jump matrix on the segment [— s, s] can be factorized using the unimodular 
fundamental solution ^(A) of (jl.lOp corresponding to the choices (|2.ip and (|2.2p 

m (I + #n(A)-02i(A) -^n(A) 
1 j " v. ^ii(A) l-^ii(A)V2i(A), 

V'll(A) "012(A)\ A -A / -022(A) -012(A)" 



021 (A) ^22(A);V° 1 / V-^lW ^ll(A) 



This motivates the first transformation of the RHP. 

3. First transformation of the RHP 
We make the following substitution in the original Y-RHP 

X(\) =Y(\)V(\), XeC\[-s,s}. (3.1) 

This leads to a RHP for the function X(X) 

• X(X) is analytic for A € C\[— s,s] 

• The following jump holds 



X+(X) = X_(A) (I A Xe[-s,s] (3.2) 
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• As A —7- ±s, we have 

A (A) = 0(ln(ATs)) 

• At infinity, 

1(A) = (/ + 0(A _1 ))tf(A), A^oo 

We will now use the Stokes phenomenon (|2,ip of ^(A) and introduce more cuts to the 
Riemann-Hilbert problem in order to uniformize its behavior at infinity. 

4. Second transformation of the RHP - uniformization 

Let 



A(A) = A(A) < 



I- 


X £ tli 


S3, 


X £ &2 


S3S4, 


xgCi 3 


S3S4SQ, 


X G CI4 



with 

Ti = {A G C : arg(A - s) = |}, T 3 = {a G C : arg(A + s) = y }, 

r 4 = {AGC: arg(A + s) = -y}, T 6 = {a G C : arg(A - s) = -^}, 
then A (A) satisfies the following RHP, depicted in Figure [2] 




Figure 2. Jump contours of the A-RHP 

• A (A) is analytic for A G C\([— s, s] U \J k T k ) 

• Along the infinite rays Tfc, A(A) has jumps described by the Stokes matrices 

A+(A) = A_(A)5 fc , A G Tfc, 
whereas on the line segment [— s, s] we have the following jump 

A + (A) = A_(A)(_°. A G [s, s]. 
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In a neighborhood of the endpoints A = ±s, 

(i, ag h h 



X(X) = X(X) 



1 -^ln^ 
1 



S3, A G 2 , ^ ^ 

S3S4, A G 
, S3S4S6, A G 0,4, 



where A (A) is analytic at A = ±s and the branch of the logarithm is fixed by 
the condition —ir < arg < 7rQ 
• As A —7- 00 the following asymptotical behavior holds 

X (A) = (/ + ^ + 0(A' 2 )) (l + ^ + 0(A" 2 )) e -(f* 3 +^V3 (4.2) 
with 

i?M 1 f-iv u\ . , 2 2 4 2 

1 2 V u iv J ' v XJ ■> x 

As we are going to see in the next sections the latter A-RHP can be solved asymp- 
totically by approximating its global solution with local model functions. Before we 
start this analysis in detail, we first connect the solution of the A-RHP to the Fredholm 
determinant det(I — K-pn). 

5. Logarithmic derivatives - connection to A-RHP 

We wish to express logarithmic derivatives of the determinant det(I — i^pn) hi terms 
of the solution of the A-RHP. To this end recall the following classical identity 



d 
ds 

In our situation 

dKp i 

ds 

where, by definition 



lndet(7 - K Pn ) = -trace ^(J - ^Pll) _1 ^^J • (5-1) 
-(A, fj,) = K Pll (X, fi) (£(// - a) + S(fx + a)) , 

s 

J 5(wt s)f(w)dw = f(±s) 



-trace ( (/ - Kpu)' 1 ^™") = _i?( S) s ) _ J?(_ S) _ s ) 



and therefore 

-t.rarp f(T— Kmr')- 1 - 

ds J 

with R(X, n) denoting the kernel (see (|2.4p ) of the resolvent R = (I — Kpu) Kpu. The 
latter derivative can be simplified using the equations (see (|2.8|) ) 

h(X) = -h 2 (X), f 2 (X) = h 1 (X) 



The local behavior (|4.ip of X(X) at the endpoints ±s can be derived directly from the a-priori 
information X(X) = 0(ln(A =p s)), A — s> ±s and the jump condition (|3.2|) . 
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as well as the identity det Y(X) = 1, which is a direct consequence of the unimodularity 
of the jump matrix G(X) and Liouville's theorem. We have, 

R , x , = F 1 (X)H 1 ( f ,) + F 2 (X)H 2 ^) = F 1 (A)F 2 ( j »)-F 2 (A)F 1 ( j ») ^ 
A — n A — fi 

Since R(X, //) is continuous along the diagonal A = fx we obtain further 

R(s,s) = F{(s)F 2 (s) - F^F^s), R(-s,-s) = F{(-s)F 2 (-s) - F^-s^-s) (5.2) 
provided Fi is analytic at A = ±s. One way to see this is a follows. Use the connection 



X(A) = Y(X)t>(X) < 



(I, AG Oi, 

03D4, A G iZ3, 
k S3S4S6, A G ^4, 

and (1221) 

F(A) = X(A)(5(A))- 1 (^(A))- 1 /(A) = X^A))- 1 -^ Q 

as well as (|4.ip to derive the following local identity 

F(A)=X(A)(J -i';fe) S (A) W )-'-i=(;)=i(A)^(;), (5.3, 

valid in a vicinity of A = ±s. But this proves analyticity of F(X) at the endpoints and 
as we shall see later on, (|5.3p is all we need to connect (|5.ip via ()5.2p to the solution of 
the X-RHP. We summarize 

Proposition 1. The logarithmic s-derivative of the Fredholm determinant ()1.12p can 
be expressed as 

-^Indetil - K PU ) = -R(s,s) - R(-s,-s), (5.4) 
ds 

R(±s, ±s) = F{(±s)F 2 (±s) - F 2 '(±s)F 1 (±s) 
and the connection to the X-RHP is established through 

where X(X) is analytic in a neighborhood of X = ±s. 

Besides the logarithmic s-derivative we also differentiate with respect to x 

d In det(J--Kpn) = -trace ((I — ifpn) -1 C ^ pn ^ 



dx \ dx J 

In our situation the kernel depends explicitly on x, since (see e.g. [T7] ) 

^ = U(\,x)% U(X,x) = -iXa 3 + i(® u JJ 
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we have 

dK <h; 1 ^' ^ = 2^(^ 21 ( A ' X )^ 11 (' U ' X ) + ^21(^,^)^11 (A, x)) 
= i(/ 2 (A)^(M)-/i(A)fei(/i)) 

and with ((231) 



S 

-trace ((J - K^T 1 ^^) = | (F 2 (X)h 2 (X) - F 1 (A)/ ll (A)) dA. 

— s 

On the other hand the Cauchy integral (|2.6|) implies 

s 

Y(A) = / + ^ + 0(A^ 2 ), A — » 00; m 1 = J F(w)h t (w)dw 

—s 

SO 

— lndet(I — Kpn) = i(m 1 — m/), mi = (m^ 7 ) 

and the connection to the AT-RHP is established via (|4.2p . Again we summarize 

Proposition 2. T/ie logarithmic x-derivative of the given Fredholm determinant can be 
expressed as 

4- bidet (J - Kpn) = i(X\ l - Xf) - v (5.5) 
dx 

with 

X(\) ~ (/+ ^1 + 0{X~ 2 ))e-^ x3+xX ^, A -> 00; Xi = (Af ). 



A 

As already indicated in the introduction, Proposition [1] and [2] are enough to determine 
the large s-asymptotics of lndet(7 — -Kpn) up to the constant term. 

6. Differential equations associated with the determinant det(7 — Kpn) 

Our consideration relies only on the underlying Riemann-Hilbert problem. Neverthe- 
less, before we go further ahead in the asymptotical analysis, we would like to take a 
short look into the differential equations associated with the AT-RHP. 

To this end we notice that the AT-RHP has unimodular constant jump matrices, thus 
the well-defined logarithmic derivatives X\X (\), X S X~ 1 (X) and X X X (X) are ratio- 
nal functions. Indeed using (|4.ip as well as (|4.2p we have 



dX 
dX = 

where 



... 1 b\ (d e \ A B 
4zAV 3 + AiX[ „ . + . , + 



X = A(X,s,x)X (6.1) 



c OJ \f -dj X-s X + s_ 

A =-i^(°o J)(*wr'i B =-^(-»»(2 i) (*(-»"' 
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and with parameters b, c, d, e, f which can be expressed in terms of the entries of mi and 
mf M (see P~2|) ). Moreover 



dX 
ds 



A 



B 



and also 



dX 

dx 



iAo"3 + i 



A + s 

b 
-c 



X = Q(\,s,x)X 



X = n(A, s, x)X. 



Hence we arrive at the Lax-system for the function X, 

A(\,s,x)X 



( dX 



^ = @(X,s,x)X, 

H = n(A, s ,x)x 

Considering the compatibility conditions of the system, 

A s -@x = [6, A], A x - n A = [n, A], 



e x -n, = [n, e] 



(6.2) 



we are lead to a system of eighteen nonlinear ordinary differential equations for the 
unknown quantities b, c, d, e, f and the entries of A and B. Since it is possible to express 
the previous derivatives of In det(7 — -Kpn) solely in terms of the unknowns b, c, d, e, f, A 
and B, one could then try to derive a differential equation for the Fredholm determinant 
(|1.12p using (|6.2|) . We shall devote to these issues our next publication. 

7. Second transformation of the RHP - rescaling and ^-function 
Along the line segment [— s, s] (see Figure [2|) we have that 

X + (A)=X_(A)(_V l ), (7.1) 

i.e. we face a permutation jump matrix. This behavior (cf. |llj ) motivates the intro- 
duction of the following g-function, 

4» 
~3 



, 1 3x 



1 ~ z, z 



oo. 



(7.2) 



g(z)=0(z) + O{z- 1 ), ${ z ) = i(\z i + ^ 2 - 



This function is analytic outside the segment [—1,1] and as z — > oo 

4 
3 

g+(z) + g_(z) = 0, ze[-l,l}. 
T{z) = X{zs)e si9{z)cT \ z e C\([-l, 1] u Jr 



Also, 
We put 



(7.3) 
(7.4) 



and, taking into account (|7.3p . are lead to the following RHP 
• T(z) is analytic for z £ C\([— 1, 1] U U fe r fc ) 
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• The jump properties of T(z) are given by the equations 



ii 



T + (z)=T„(z) 



-i 
-i 



In a neighborhood of the endpoints z = ±1 



T( z )e- si 3^ = X{zs) (I 



Z7T Z + l 



r iii 




, zer k . 




I- 


A G Qi 


S3, 




S3S4, 


A G A 3 


S3S4SQ, 


AgQ 4 



(7.5) 



• As z — > 00, we have T{z) = I + 0(z 1 ) 

Let us analyse the behavior of the jumps along the infinite branches 1^ as s — > 00. To 
this end consider the sign-diagram of the function Re g{z), depicted in Figure where 
x is chosen from a compact subset of the real line and s > is sufficiently large. 




Figure 3. Sign-diagram for the function Re g{z). Along the solid lines 

7T I 27T 

■3' 31 3 



Re g(z) = and the dashed lines resemble arg z = ±5, ±% 



Since Re g(z) is negative resp. positive along the rays Ti,r3 resp. T4,Tq, 



uniformly on any compact subset of the set (|1.15|) and the stated convergence is in fact 
exponentially fast. From (|7.6|) we expect, and this will be justified rigorously, that as 
s — > 00, T(z) converges to a solution of the model RHP, in which we only have to deal 
with the constant jump matrix on the line segment [—1,1]. Let us now consider this 
model RHP. 
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1 \ !/ 4 
Z + 1 N 1 



8. The model RHP 

Find the piecewise analytic 2x2 matrix valued function M{z) such that 

• M(z) is analytic for z G C\[— 1, 1] 

• Along [—1,1] the following jump condition holds 

M+(z) = M-(z){? i - i ), «e[-l,l] 

• M(z) has at most logarithmic singularities at the endpoints z = ±1 

• M(z) = I + 0(z~ 1 ), z^oo 

A solution to this problem can be obtained explicitly via diagonalization (cf. [11] ) 
with 

^ i . ! 

and (f^j) 1//4 is defined on C\[— 1, 1] with its branch fixed by the condition (f^j) 1//4 — > 1 
as z — > oo. 

9. CONSTRUCTION OF A PARAMETRIX AT THE EDGE POINT Z = +1 

In this section we construct a parametrix for the right endpoint z = +1. We start 
from the local expansion 

= ^p* (f + 0) V^T(l + 0((z - 1)V2)) , *->!,-*< arg(z - 1) < vr 

and the singular endpoint behavior (|4.ip 

X(z) = 0(]n(z-1)), z^l. 

Both observations suggest (cf. [11]) to use the Bessel functions H^^) and Hq 2 \^) for 
our construction. This idea can indeed be justified rigorously as follows. First we recall 
that the Hankel functions of first and second kind are unique independent solutions to 
Bessel's equation 

zw" + vJ + w = 

satisfying the following asymptotic conditions as £ — > co and — ir < arg £ < ir (see [T]) 

Secondly H^\^),Hq(Q satisfy monodromy relations, valid on the entire universal 
covering of the punctured plane 

H^iCe™) =-i4 2) (C), ^ 2) (Ce") =^ 1) (C) + 2^ 2) (C), ^ 2) (Ce"") = -^(C) 

(9.1) 
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and finally the following expansions at the origin are valid (compare to (|4.ip ) 

(C) = «o + 01 In C + a 2 C 2 + a 3 C 2 ln( + (C 4 In C) , C -> (9.2) 
with coefficients given as 

2ry E 2z 2i i 1 i i 

o = H In 2, ai = — , a 2 = — (1 -7s) - - + — ln2, a 3 = - — 

7T 7T 7T Z7T 4 Z7T 27T 

where is Euler's constant and the expansion for Hq (£) is up to the replacement 
dj i— > a, identical to (|9.2p . Remembering all these properties we now introduce the 
following matrix-valued function on the punctured plane £ £ C\{0} 



Pbe(C) = 



r(2) 



r(!) ' 



e 



-7T < arg < 7r. (9.3) 



Using the behavior of Hq 1 ^ (£) and i?Q 2 ^ (£) at infinity we deduce 



Pbe(C) 



15* 



1 1 
-1 1 

6 



1 + 



8VC V 2 1 



2\ 3 
+ 



128C 



1 -4 
-4 1 



x +0(C 



1024C 3 / 2 V~ 6 

as C - > oo and — ir < arg £ < ir. Let us now assemble the following model function 



pffi (0 = { 



Pbe«) 



Pbe(C) 



1 

-% 1 
1 % 

1 



which solves the RHP depicted in Figure [H 



arg C G (f ,tt), 

arg C G (-tt,-|; 
arg C €(-£,§). 



(9.4) 



1 o - 

-! 1, 



5i 



-i" 
-i , 



(o 1*)- 5 6 

Figure 4. The model RHP near z = +1 which can be solved explicitly 
using Hankel functions 

More precisely, the function P§§ (0 possesses the following analytic properties. 
• (0 is analytic for C G C\{arg £ = -tt, -§, §} 
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• The following jumps hold 

W(C)) + = W(C))_(J 7), arg C = -^ 

And for the jump on the line arg £ = ir we notice that the monodromy relations 
imply 

H?\Vt + ) = ^ 2) (V / C_e-)=< ) (7C_)+2^ 2) (7C_) 

(^ 2) )'(Vc+) = e-(<))'(7c_) + 2e-(^)'(yc_) 

and 

^(VC*) = ^ 1) (>/C-e" < )=-flj a) (>/C-) 

« } )'(^ + ) = (^ 2) )'(VU- 

Therefore 



and hence 



^#f(C)) + =(P#f(C))_(_°. "'), argC = vr, 



In order to determine the behavior of -Pgjf (C) at infinity we make the following 
observations. First let arg £ S (f i 71 ") an d consider 



Observe that Refz-y^C) < 0, hence the given product approaches the identity 
exponentially fast as £ — )• 00. Secondly for arg £ E (— tt, — 5) 

and in this situation Re( — iy/C) < 0, so again the product approaches the 
identity exponentially fast as ( — > 00. Both cases together with the previously 
stated asymptotics for Pbe(C) imply therefore 



15i / 1 6 



1024C 3 / 2 

as C —7- 00 in a whole neighborhood of infinity. 



e^ CT3 , (9.5) 
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The model function Pg§ (C) will now be used to construct the parametrix to the solution 
of the original X-RHP in a neighborhood of z = +1. We proceed in two steps. First 
define 

((z) = -s 6 g 2 (z), \z - 1| < r, -vr < arg ( < vr (9.6) 

or respectively 

This change of variables is indeed locally conformal, since 

«*) = ^ (§ + |j) V " 1)(1 + 0(* - 1)), k-l|<r 
and it enables us to define the right parametrix U(z) near z = +1 by the formula: 

U(z) = Br(*)y ^1 e-**J$f (CWK 3 ^ 3 , k - II < r (9.7) 
with C(z) as in (|9.6p and the matrix multiplier 

By construction, in particular since B r {z) is analytic in a neighborhood of 2 = +1, the 
parametrix U(z) has jumps along the curves depicted in Figure [5j and we can always 
locally match the latter curves with the jump curves of the original RHP. Also these 
jumps are described by the same Stokes matrices as in the original T-RHP. Furthermore, 
and we will elaborate this in full detail very soon, the singular behavior of U (z) at the 
endpoint z = +1 matches the singular behavior of T(z): 

U{z) = 0(\n(z- 1)), |z-l|<r. (9.9) 

Hence the ratio of T(z) with U(z) is locally analytic, i.e. 

T(z) = N r (z)U(z), \z - 1| < r < -. (9.10) 




Figure 5. Transformation of parametrix jumps to original jumps 
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Let us explain the role of the left multiplier B r (z) in the definition ([9 
that 
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Observe 



B r (z)((z) 



-ct 3 /4. 



1 1 1 



2 \l -1 



M(z). 



This relation together with the asymptotic equation (|9.5p implies that 



U(z) 



+ 



1 1 

15i 



1 

1 -1 



/ + 



-1 



+- 



1024C 3 / 
/ + 



-6 -l) + ^~ 



8^C 



1 A 1 



-2\ _3_ 
1 ) + 128C 



2 U 



1 1 
1 -1 



1 

M(z) 



128C 



■(5/? 2 -7/3" 2 ) 5/3 2 + 7/T 2 \ , 2 v 



2048C 3 / 2 V"( 5 / ? + 7 ^ ) W ~ 7 P~ 



M(z) 



(9.11) 



as s — > 00 and < n < \z — 1| < r 2 < 1 (so |£| — > 00). Since the function £(2) is of order 
0(s 6 ) on the latter annulus and f3(z) is bounded, equation (|9.11|) yields the matching 
relation between the model functions U(z) and M(z), 

U(z) = (I + o(l))M(z), s^oo, < n < |z - 1| < r 2 < 1 

which is crucial for the successful implementation of the nonlinear steepest descent 
method as we shall see after the next subsection. This is the reason for chosing the 
left multiplier B r (z) in (p^) in the form (f9T8j) . 



10. Construction of a parametrix at the edge point z 



-1 



The model RHP near the other endpoint z = — 1 can be introduced in a similar way 
as we did it in the last section. First we consider on the punctured plane £ € C\{0} 




-e 2 



< arg C < 2vr 



1^3/4 



-1 -1 

-i i 



1 + 



8VC 



+- 



15 



-1 6 



; +o(c- 



-1 2 
-2 1 



+ 



128( 



-1 
-4 



1024C 3 / 2 V~ 6 

as C — > 00 and < arg ( < 2tt. Next, instead of (|9.4p . define 



pRH 



(0 



Pbe(Q 

Pbe(C) 
Pbe(C), 



1 

-i 1 

1 i 
1 



arg C G (0, f ), 



argCG (^,2vr), 
arg Ce (¥,¥)• 



(10.1) 
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which solves the model RHP of Figure [6) More precisely, the function Pgjf (C) has the 
following analytic properties 




Figure 6. The model RHP near z = — 1 which can be solved explicitly 
using Hankel functions 

. P§§ (C) is analytic for C € C\{arg C = f , 7 f , 2tt} 

• We have the following jumps on the contour depicted in Figure [6] 

(Pff(C)) + = (Pff(C))_Q J), arg C =f 

(P#f(C)) + = (P#f(C))_(j j), arg C =^ 

and on the line segment arg £ = 27r 

jtfV'VCj = ^ 1) (e-^ v ^_e-)=^(e-iv / C-)+2^ 1) (e-^7C_) 
(^V 1 ^) = e-(^)'(e-f7cJ+2e-(^ 1) ) / (e"^7Cj 



as well as 
hence 

(^f(C)) + = (^(0). 



arg £ = 27r. 



A similar argument as given in the construction of Pgjf (C) implies 



^(0 



0-3/4 



/ + 



+- 



15 



-1 6 
-6 1 



8VC 

+ o(r 2 ) 



-1 2 
-2 1 



+ 



128C 



-1 

-4 



(10.2) 



1024C 3 / 2 

as C — 7- 00, valid in a full neighborhood of infinity. 

Again we use the model function Pg§ (C) in the construction of the parametrix to the 
solution of the original X-RHP near z = —1. Instead of (|9.6p 

(( z ) = s 6 g 2 (z), \z + 1| < r, < arg C < 2vr (10.3) 
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or equivalently 

This change of the independent variable is locally conformal 



32s 6 /3 3x x 2 



\ (z + l)(l + 0(z + l)), \z + l\ < 



9 V 2 4s2 , 

and allows us to define the left parametrix V{z) near z = — 1 by the formula: 

V(z)=B l (z) 1 -^ ^y|p|f(C(,)) e ^V3, |* + l|< r (10.4) 
with the matrix multiplier 

(10.5) 

Similar to the situation in the last section, V{z) has jumps on the contour depicted in 
Figure [7] which are described by the same Stokes matrices as in the original T-RHP. 

e ~s 3 g(z)a 3 g^ e s a g(z)cr 3 

(- < 




3 S 4 e 8 

Figure 7. Transformation of parametrix jumps to original jumps 

Also here, as we shall see in detail in the section [131 the singular behavior at z = — 1 
matches: 

V(z) = 0(hx(z+ 1)), |« + l|<r (10.6) 
Hence the ratio of parametrix V(z) with T(z) is locally analytic 

T(z) = N l {z)V{z), |z + l|<r<^ 
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and the left multiplier (|10.5j) in (|10.4p provides us with the following asymptotic matchup 
between V{z) and M{z): 



V{z) 



+ 



1 1 
1 -1 

15 



P(z) 



1 (\ 1 



0-3 

2 I 
1 6 N 



/ + 



1 



-1 2 



+ 



o(r 2 ) 



1 + 



+ 



1024C 3 / 2 V" 6 1 

1 / /3" 2 -3/3 2 
16VC *v-(r 2 + 3/3 2 ) 
15 / 5/T 2 -7/3 2 
(5/T 2 + 7/3 2 ) 



8VCV- 2 1 
1/1 1 
2 ^1 -1 

/3" 2 + 3/3 2 
(/3- 2 - 3/3 2 , y 

5/T 2 + 7/3 2 
-(5/T 2 -7/3 2 ) 



+ 



128C 
1 1 
1 -1 

1 - 
4 - 



M(z) 



128C 

+ o(r 2 



Af(z) (10.7) 



as s 



2048C 3 / 2 

00 and < 7*1 < |£ + 1| < r2 < 1, thus 

V{z) = (I + o(l))M(z), s^oo, 0<n< |* + 1| <r 2 <l. 

At this point we can use the model functions M(z), U{z) and V{z) to employ the final 
transformation. 

11. Third and final transformation of the RHP 
In this final transformation we put 



R(z) = T(z) < 



' (U(z))-\ 
(Viz))' 1 , 
{ (M(z))-\ 



\z — 1| < r, 
\z + l\ < r, 
\z =p 1| > r 



(11.1) 



where < r < -7 remains fixed. With C T and C\ denoting the clockwise oriented circles 
shown in Figure [HJ the ratio-function R{z) solves the following RHP 



73 




74 




Figure 8. The jump graph for the ratio-function R{z) 



• R(z) is anlytic for z G C\{C r U Q U \J k 7/J 

• Along the infinite branches 7 fc 



R + (z) = R^z)M(z)e~ s3 ^S k e sS ^{M(z))-\ z G 7 fc, 
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and on the clockwise oriented circles C r \ 

R + (z) = R-(z)U(z)(M(z)y 1 , zeC r 
R+(z) = R.(z)V(z)(M(z)y\ zed 

• R(z) is analytic at z = ±1. This observation follows from (|9.9p and (|10.6p . which 
will be proved in section [T3j 

• In a neigborhood of infinity, we have R(z) — > I 

We emphasize that, by construction, the function R(z) has no jumps inside of C r and C\ 
and across the line segment in between them. In order to apply the Deift-Zhou nonlinear 
steepest descent method for the ratio- RHP, all its jump matrices have to be close to the 
unit matrix, as s — > oo, see [H]. Due to the triangularity of all Stokes matrices Sk, the 
jump matrices corresponding to the infinite parts (J fc 7^ of the -R-jump contour are in 
fact exponentially close to the unit matrix 

||Me-^(>^ fe e^> 3 (M)" 1 - I\\ L ^L°°{ lk ) < cie"^ 2 ^ (11.2) 
emanating from C r j as s — > 00 with constants Cj > whose values are not important. 
Moreover, by virtue of (|9.1ip . f7(z) (M(2i)) 1 approaches the unit matrix as s — > 00 

WU(M)- 1 - I\\ L 2 nL oo {Cr) < c 3 s~ 3 (11.3) 

and from (fTtTTlh also V(z)(M(z)y 1 

Hy(M)" 1 - I\\ L 2nL°°( Cl ) < c 4 s" 3 , s -> 00. (11.4) 

All together, with Gr denoting the jump matrix in the latter ratio-RHP and the 
underlying contour 

\\ G R ~ m^nL^iEn) < cs" 3 , (11.5) 
uniformly on any compact subset of the set (|1.15p . The latter estimation enables us to 
solve the ratio-RHP iteratively. 

12. Solution of the RHP for R(z) via iteration 

The RHP for the function R(z), 

• R(z) is analytic for z £ C\S^j 

• Along the contour depicted in Figure [8] 

R+(z) = R-(z)G R (z), z£S fi 

• At infinity, R(z) = I + 0(z~ 1 ),z 00 
is equivalent to the singular integral equation 

R_( Z )=I +7 L f R_ (y;) (G R (w) - I) — (12.1) 

ZTTl J W — Z- 

and by standard arguments (see [14] ) we know that for sufficiently large s the relevant 
integral operator is contracting and equation (|12.1|) can be solved iteratively in L 2 (X#). 
Moreover, its unique solution satisfies 

\\R- ~ I\\l*(x r ) < cs~ 3 , s^oo. (12.2) 
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It is now time to connect the latter information with (|5.4p and f)5 .5 j) to determine the 
large s asymptotics of det(i — K-pu) up to the constant term. 

13. Asymptotics of lndet(7 — -Kpn) - preliminary steps 

Let us recall the transformations, which have been used in the asymptotical solution 
of the original Y-RHP 

Y{\) i y 1(A) ^ X(X) ^ T(z) i— > R(z). 

Thus, in order to determine lndet(I — -Kpn) via Proposition Q] we need to connect X(±s) 
and X'(zLs) to the values i?(±l) and i?'(±l) of the ratio-function. This can be done as 
follows: From and (|7.5|) for \z — 1| < r 



R(z)U(z)e- s3 ^ z)r73 = X{zs) 
and for | z + 1 1 < r 

R{z)V{z)e~ s39 ^ as = X{zs) 



2ir z+1 



1 7 



/, AG Hi, 

53^4, A e 



-S3, A e r^2, 
S354, A e S^3- 



(13.1) 



(13.2) 



This shows that the required values of X(±s) and X'(±s) can be determined via com- 
parison in either (|13.ip or (|13.2|) . once we know the local expansions of U(z), respectively 
V(z) at z = ±1. Our starting point is (|9.2f) 



pRH 



(0 



ao + flnC i(a + ^lnC) N 
a\ ai 

9 -"l + -V^ r / (a i + ! ln 9,)+0(C 2 lnC), (13.3) 
2a 2 + a 3 + a 3 In C 2a 2 + a 3 + a 3 In C / 



as C — )• and — § < arg £ < ^. The latter expansion together with the changes of 
variables C = C( z ) = ~~ s 6 5 2 (^) an d A = zs implies for — | < arg (A — s) < ^ 

Pbe{((z)) =n(ln(A-s)) +(A-s)n(ln(A-s)) + 0((A - s) 2 ln(A - s)) , A^s 

with the matrix functions = (fiy) and II = (Ily) being determined from (|13.3p . Now 
we combine the latter expansion with (|9.7p and (|11.1|) 

i?(z)^(2)e- s3 ^|^ A=J R(z)i? r (z)^J|e-fP##(C(z))| 2=A 

s » s 

= fl(l)B r (l)^y|e- < in(MA - «)) + (A - a) j/2(l)fl r (l) 

xf yfe^fnlMA-,)) + (R>(l)B r (l) + R(l)B> r (l))°f ^^ "(MA- ^)) 
+0((A-s) 2 ln(A-s)), A^s, -| < arg (A - s) < | 
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and similar identities hold for — ir < arg (A — s) < — ^ and ^ < arg (A — s) < ir, they 
differ from ()13.3|) only by right multiplication with a Stokes matrix (see (|9.4p ). On the 
other hand the right hand side in ()13.ip implies for — 5 < arg (A — s) < ^ 



X(s) + (A - s)X'(s) + 0((A - s) 2 ) 



2tt 111 A+s 



1 

i 1 



/iu(s) ^ln(2 S )A > 11 (s)+X 12 ( S ) 
VX 21 (s) ±ln(2s)X 21 (s)+X 22 {s) 



+(X-s) 



X' u (s) ±(ln(2s)X' u (s) + i- s Xii(s)) + X[ 2 (s) 
X' 21 (s) i F (ln(2 S )X 21 ( S ) + ^X 21 ( S )) + X 22 ( S ) 



i-SFM*"*) -^ ln ( A - s ))+0(( A - s ) 2 ln(A- S )), A 



and thus by comparison of the left hand side and right hand side in (|13.ip 
X n (s) = ^e-^(R(l)B r (l)) lv X 21 (s) = ^-e~^ (fl(l)B P (l)) 



21' 



(13.4) 



Although (|13.4j) was derived from a comparison in the sector — f < arg (A — s) < -|, the 
multiplication of f 1 1 3 . 3 j) with the right Stokes matrix in the other sectors as well as the 
use of appropriate Stokes matrices in (|4.ip show that (|13.4p follows from comparison in 
a full neighborhood of A = +s. Comparing terms of 0((A — s) ln(A — s)j we also derive 



X' u (s) 



is 5 / 3 3x\ 2 



T\ 2 + 4^ j ((W^)ll + 2i W) fi r(l)) 12 



--(^(l)B r (l) + 12(1)4(1)) 



ii 



(13.5) 



and 



*21 (*) 



e 4 



sV3 3x 



9 V2 



+ 



—A [(R(l)B r (l)) 21 + 2i(R(l)B r (l)) 22 



--{R'(l)B r (l) + R(l)B' r (l)) 21 
All together we have the identities 



(13.6) 



Fi(s) 



Xu(s), F 2 (s) 



X 2 i(s), F{(s) 



X' 21 (s) 



related to the solution of the ratio- RHP via (j 13.4[) . f|13. 5j) and (|13.6p . A completely 
similar analysis for the left endpoint A = — s provides us with 

Fi(s) = J(#(-l)Bz(-l)) 12! F 2 (-s) = l -(R(-l) Bl (-l)) 22 (13.7) 
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and 



+-(i?(-i)^(-i)+i?'(-i)^(-i)) 



12 



(13.8) 



as well as 

We can now derive 



8s V 3 3x 



9 V2 



+ £2 ) ((^(- 1 )^(- 1 )) 2 2 + 2(i?(-l)^(-l)) 21 



(i?(-l)^(-l) + J R'(-l)^(-l)) 



22 



(13.9) 



R(s,s) = F{(s)F 2 (s) - F^F^s) 



4s V 3 3x 
2 + 4s2 



(J?(l)J3 r (l)) 11 (fl(l)5 r (l)) 2a - (i2(l)J3 r (l)) 21 (fl(l)J3 r (l)) 



12 



+ 



4s 



(^'(1)^(1) +i2(l)S;(l)) n (i2(l)5 r (l)) 21 - {R'(l)B r (l) + R(1)B>(1)) 



21 



(i?(l)S r (l)) 



11 



as well as 



R(-s, -s) = F{(-s)F 2 (-s) - F^-s^i-s 
4s 5 ^3 3s y 



9 V2 4s 2 
x(fl(-l)B,(-l)) 12 



(22(-l)B,(-l)) 11 (i2(-l)B l (-l)) 22 " (^(-1)^(-1)) 21 
l r, 



4s 



(i?(-l)^(-l)+i?'(-l)A(-l)) 



12 



(i?(-l)^(-l)) 22 -(i?(-l)^(-l) + J R'(-l)^(-l)) 22 (i?(-l)^(-l)) 



12 



We make the following observation 

Proposition 3. R{z) is unimodular for any i£M, i.e. det-R(z) = 1. 

Proof. From ([93D we obtain that detP^f (C) = ~, hence det [/(z) = 1. Similarly 
det -Pgjf (C) = ~~ — leading to det V(z) = 1. Thus the ratio- RHP has a unimodular jump 
matrix Gr(z) which implies by normalization at infinity deti?(z) = 1. □ 

Applying the latter Proposition, one checks now readily 

(i2(l)B P (l)) 11 (i2(l)S r (l)) 22 - {R(l)B r (l)) 21 (R(l)B r (l)) 12 = -2 

and 

(i2(-l)J3 I (-l)) 11 (i2(-l)fl J (-l)) 22 " (^(-l)A(-l)) 21 (i?(-l)^(-l)) 12 = "2. 
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We combine these two identities with the values of B' r (l) and B'A— 1) to deduce 



R(s,s) 

as well as 
R(-s,-s) 



2is v 



3x \ 



8s 5 /3 3x\ 
-(R 21 (l) + R 22 (l))(R n (l)+Ri2(l)) 



{R' u (l) + R' 12 (l)){R 21 (l) + R 22 (l)) 



8s 5 /3 3x x 



2is 2 /3 3x 
T 2 + 4? 



(i4(-i)-# 12 (-i)) 



(13.10) 



(13.11) 



(i?2l(-l) " i?22(-l)) " (i^i(-l) " # 22 (-l)) (flll(-l) " i?12(-l)) 



14. ASYMPTOTICS OF lndet(7 - Kpn) UP TO THE CONSTANT TERM 



The stated asymptotics (I1.14p without the constant term is a direct consequence of 
Proposition [1] and [2j We first need to compute i?(±l). From the integral representation 
and estimates (fTT31) . (fl~!T2"]l 

= / + i / (G R ( W )-l)^- 1 +0(s-*)=I + 0(s-3), 



OG 



so 



i?(s,s) 



8s 5 /3 3x 



2is 2 /3 3.x 
9 \2 + I^) + ^T\2 + I^ 

+R' 12 (1) - R' 21 (l)] +0(s" 4 ) 



R' n (l) - R! 22 (l) 



and 



8s 5 /3 3xV 2is 2 /3 3x 



i?ii(-l)-i? 22 (-l) 



- J R / 12 (-1)+ J R' 21 (-1)] +0( S " 4 ), s^oo. 
In order to compute the values R'(±l) one uses (|9.11[) and (jlO.71) 



R>(±1) 



1 
1 



- /) 



dw 



(w^l) 2 



dw 



/3 2 - 3/3- 2 -(/3 2 + 3/3- 2 ) 
2ttz7 l 6 ycR V/ 3 ^ 3 / 3 " 2 -(/3 2 -3/3" 2 ); (»Tl) 2 



+- 



2vrz 



^-2 _ 3/3 2 r 2 + 3/3 2 

-2 i oo2\ ( o-2 OQ2-- 



i6Vch + m -{p- 2 - mj (w t i) 



dw 



ASYMPTOTICS OF A PII FREDHOLM DETERMINANT 28 

with the local variables given in (|9.6p . (|10.3p : 

(3 2 {w) 3/ 2 1 3x\~ l 1 0~ 2 (w) 3/ 2 1 3x\ _1 1 
w € C r : —7=== = — r ( to +— + — ^ I 7, , . . = -r-Q I «r+^ + 



/3 2 H 


3 




~ 4s 3 




Si 




~ 4s 3 



2 4s 2 y to-i' ycH 4s3 v 2 Is 2 ) w + 1 

9 1 3a; \ 1 1 /3 _2 (u>) 3i / , 1 3x\ _1 1 



In the end residue theorem leads us to 

'3 3x \ 3i ( — 1 =pl 



and 



256s 3 1 

+ 64^2 + 4s 2 J Uf + 4) f+lfc 
3i /3 3x V 2 /'-l ±l\ /-,/ -e\ 



8o5 / q q„, \ 2 q 

i?( s , s ) = i?(- s ,- s ) = ^(j + 0j +± + ( S - 3 ), s^oc. (14.1) 

Combining (|14.ip with (|5.4p we have thus derived the following asymptotics 
d 3 

— lndet(/-Jfpn) = -4s 5 -4xs 3 -x 2 s- — + 0(s" 3 ), s -> 00 (14.2) 
as As 

and the error term is uniform on any compact subset of the set (|1.15|) . 

Opposed to (|14.2p we are going to determine det(7 — K-pu) now via Proposition [2] 

4- bidet (J - Kpu) = i(X\ l - Xf) - v 
ax 



where 

To this end remember the definition of /3(z) and g(z), hence as z — > 00 



Xf = Um (\(X{\y^ x +xX > as - 7) ..J , i = 1, 2. 



mm =z+i(; ;)+o(,- 2 ), e -w- fl wH =/ + ^(i+ j)a 3 +o(,- 2 ), 



which gives 

is 3 rc \ 1/^0 1^\ i 
"2 



already neglecting exponentially small contributions in the last equality. The latter 
integral can be evaluated in a similar way as we did it during the computation of (|14.2p . 
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we end up with 



i 

2^ 



C r ,l 

i 

2^ 



R.(w){G R { W )-l))^dw = 
R-(w){G R (w) - 1)} dw, 



3i I 3 3x 
+ 



32s d V 2 4s^ 



i.e. together 
d 



dx 



lndet(I - Kpii) = 2is 



it 



1 + 4H 



3i f 3 3x 
32t \2 + 4? 



+ 0( S " 5 ) 



1 



s — oo, 



(14.3) 



again uniformly on any compact subset of the set (|1.15p . Given the asymptotic expan- 
sions (|14.2p and f 1 1 4 . 3 [) we can now determine the large s-asymptotics of lndet(I — i^pn) 
via integration 



lndet(J- K s 



s 2 -jjlns + j (y-x)u 2 {y)dy + uj + 0{s~ l ), (14.4) 



recalling that u(x) ~ Ai(x) as x — > +oo. As we see (|14.4p matches (|1.14p up to a 
general, universal constant term oo. We shall now determine this constant term using an 
approximation argument for the given kernel (jl.lip . 



15. Kernel approximation - Kpn *->■ i^csm 

Within the asymptotical analysis of the X-RHP in the past sections, one of the first 
transformations of the RHP was the (^-function transformation: it allowed us to trans- 
form the jumps on the infinite branches to exponentially small contributions and to 
solve the model problem on the line segment [—1, 1]. Then and there it was crucial that 
for x chosen from a compact subset of the real line and s sufficiently large, one always has 
that Re g(z) is negative on the infinite parts Ti,r3 in the upper halfplane and positive 
on the infinite contours T4, Tq in the lower halfplane. This fact however also holds in 
the limit x — > +00, on the other hand it fails for x — > —00: Let 

, . fl 3x~ 
Z± = ± l2 + I^ 

denote the two vertices of the curves depicted in Figure below. In case x,s > 0, they 
are purely imaginary and bounded away from zero, hence the statement on the sign of 
Re g(z) on r, follows. 



Our approach henceforth will be to study the large positive x-limit of (jl.lip . i.e. the 
large positive x-limit of the associated function ^(X,x). We begin with the following 
Riemann-Hilbert problem depicted in Figure [TOl compare [T7] 
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Figure 9. Sign diagram for the function Re g(z). In the left picture 
we indicate the location of z± as x > and in the right picture for a 
particular choice of x < 0. Along the solid lines Re g{z) = and the 
dashed lines resemble arg z = db^db^ 



S3 Si 




Figure 10. The RHP jump graph associated with the Hastings-McLeod transcendent 

• V I' 00 (A) is analytic for A € C\(|J fe i2fc) where Rk denote the rays 

R k = {\£C\ arg A = J + - 1)}, fc = 1,3,4,6 

u 6 

• On the rays Rk, the boundary values of the the function *$>°° satisfy the jump 
relation 

tf~(A) = y°°(\)S k , X £ Rk, k = 1, 3, 4, 6 

• At A = 00 the following asymptotic behavior takes place 

^ 00 (A)e i( t A3+xA),T3 =1 + 0(A- X ) 

which is connected to the given ^-function of p. lip by 

*(A,x) = *°°(A,x)5i. 

As we see, determining the large positive x behavior of ^f(X,x) therefore reduces to an 
analysis of the oscillatory ^^-RHP. However the latter RHP is very well known since 
it is used to determine the large x-asymptotics of the Hastings-McLeod solution of the 
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second Painleve transcendent given in the introduction (cf. [H]). We have in fact for 
A £ (s,s) 

^°°(X,x)e i ^ x3+xX ^-I = o( X ~ 1 / /4e ~' X ' V x^+oo 

V V4A 2 + x / 

hence 

- _1/4 _2~3/2 



(A, z) = e-^^+O^ 
^21 (A, x) = -ie^i x3+ ^+0^ 



V4A 2 + x 

1/4 -2-3/2 



V4A 2 + x 

as x —7- +oo and A £ (— s, s). Going back to (jl.lip we obtain 

/ / x^e-l^ 2 \\ 
K Pn (\,ij,) = K csin (X,n)\ 1 + , x -> +oo, A,/i £ (-s,s) 

y V V(4A 2 + x)(4,u 2 + x)J J 

(15.1) 

where 

- , , sin(|(A 3 -/i 3 )+x(A-/i)) 

tf csin A, n) = ^ ,\ / — • ( 15 - 2 ) 

Vr(A - /x) 

The latter integral kernel is a cubic generalization of the well known sine kernel 

sinx(A — fx) 
vr(A - fi) 

acting on L 2 ((— s, s); <iA). In order to compute the constant term in (|1.14p we will 
introduce a parameter t £ [0, 1] and pass from (|15.2p to 

- , s , s sin (|t(A 3 - n 3 ) + x(A - fi)) , 

^csin(A, /X M- ^csin(A, y) = ^ ^ ' / ^ (15.3) 

Our strategy is as follows. First we are going to find the large s-asymptotics of 

^lndet(/-K csin ) (15.4) 

using the Riemann-Hiibert approach of section 2. This analysis can be done indepen- 
dently of the original kernel (jl.lip , although in some details it is similar to the presented 
one. Afterwards, using uniformity of the asymptotic expansion with respect to t £ [0, 1] 
we shall integrate 
l 

— lndet(J - -FT C sin) dt = lndet(J - K csin ) - lndet(7 - K sin ); 



J dt 

o 



but since the asymptotic expansion of the sine kernel as s — > oo is known including the 
constant term, we know the large s-asymptotics of 

det(J - K csin ) 

also up to order 0(s^ 1 ), in fact 

In det(I - K csin ) = A(s,x) +lj + O(s~ 1 ), s -)• oo (15.5) 
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uniformly on any compact subset of the set (jl,15p with 

A(s,x) = -h 6 - s"x - l(sx) 2 -hns, u = -|ln2 + 3C'(-1), 
which is the statement of Theorem [2j On the other hand from (|14.4p 

oo 

In det(7 - i^pn) = A(s,x) + J (y - x)u 2 {y)dy + uj + 0(s _1 ), s ->■ oo (15.6) 

X 

hence considering (|15.1|) as well as 

oo 

lim f(y- x)u 2 {y) = (15.7) 

x— >oo J 

X 

we might conjecture that uj = ujq, and this conclusion can be justified as follows. Notice 
the following identity for trace class operators A, B (see |27j ) 



det(I — A){I — B) = det(/ - A) det (7 - B) 
which gives in our situation 

det(7-£r Pn )-det(7-7? csin ) = -det(7-7T csin )[l-det (i- (J-^ CS m) _1 (^Pii-^c S in) 
provided 

(/ - i^csin)- 1 = 7 + R csin (15.8) 

exists as a bounded operator. The latter statement will follow from the Riemann-Hilbert 
analysis given in section [23l Since from (|15.ip 

s 

(K PU f)(X) = J Km&rifinW = (i<W)(A) + (Ef)(X) 

—s 

where the trace class operator E has a kernel satisfying 

E(X,fi) = o(x~ 1/4 e-^ x3/2 ^, x^oo, (A,//) 6 [s,s] x [s,s] (15.9) 

we obtain 

det(/ - K PU ) - det (7 - K csin ) = - det (7 - 7? csin ) [l - det (7 - (7 + R csin )E^ 
and therefore 

det (7 - 7^ csin ) V / 

Now from the boundedness of 7 + 7^ cs i n as well as (|15.9|) we see that the convolution 
kernel of the operator 

(7 + R cs i n )E 
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approaches zero exponentially fast as x — > oo, thus via Hadamard's inequality 

det (/-(/ + Acsin)^) 

oo s s 

= 1 + ^2 " T / ' " / det [(^ + ^csm) E ] (Xi,Xj)dxx ■ ■ ■ dx n = 1 + O s (l), X ^ OO 
or similarly 

lndet(I--Kpii) = lndet(7 - X csin ) + o s (l), x oo. (15.10) 
We combine (|15.5|) . (|15.6p . (|15.10p and obtain from Schwarz inequality 

oo 

. I (3(s) f . o/ \ r 

wo - w < - H h / {y - x)u (y)dy 

s x J 

x 

for all x > xq and s > so, with a universal constant a and a positive function /3 = /3(s). 
Recalling (|15.7|) we first take the limit x — > oo and afterwards s — > oo to conclude 
uo = w. 

16. RlEMANN-HlLBERT PROBLEM FOR det(7 — i^csin) 

Before we start proving (|15.5|) , let us locate (|15.3|) within the framework of integrable 
Fredholm operators and derive the connection of (|15.4p with the solution of the under- 
lying Riemann-Hilbert problem. 

The given kernel (j!5.3[) is of integrable type with 

_ d\\)e^) 1 / e ^A3+,A)x x , e - i{ ltxs +xX) 

hence Lemma 2 implies the following O-RHP 

• @(A) is analytic for A £ C\[— s,s] 

• On the line segment [— s, s] oriented from left to right, the following jump holds 



• 0(A) has at most logarithmic endpoint singularities at A = ±s 

9(A) = 0(ln(A=F»)), A -> ±s 

• 0(A) — > I as A — > oo. 
Also here we can factorize the jump matrix 



^ W ~ I _ e -2i$t\*+xX) 2 ]- e3 1-1 2' 

and employ a first transformation. 
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17. First transformation of the 6-RHP 

Introduce 

$(A) = G(A)e i( t a3+xA)<T3 , AGC\[-s,s] (17.1) 
and obtain the following RHP 

• 3>(A) is analytic for A G C\[— s, s] 

• On the segment [— s,s] 

*+(A) = *-(A)(_° 1 J)' 

• ^(A) has at most logarithmic endpoint singularities at A = ±s 

• We have the following asymptotics at infinity 

$(A)~ (7 + 0(A- 1 )y(i' A3+xA ) (T3 , A ^00. 

As we are going to see in the next sections the latter RHP admits direct asymptotical 
analysis via the nonlinear steepest descent method. This analysis shows similarities to 
the analysis of the sine kernel determinant analysis (see [11]) and the analysis presented 
in the past sections of the current paper. However one major difference to (|1.11|) is 
the absence of infinite jump contours in the given $-RHP, hence we should not start 
our analysis from the X-RHP in section 2 and use the previously discussed large x- 
approximation ^f°°(X,x). 

Before we start the asymptotical analysis, let us first express (|15.4p in terms of 6(A). 



18. THE LOGARITHMIC t-DERIVATIVE of det(J - K csin ) 

The following derivation is a special situation of the general case given in [23]. We 
start from (|5,ip 

s 

|lndet(/-ir csin ) = -J f(l-K c8in y 1 ^-y\,X)dX 

—s 
s 

= -f ((i +Rcsin )^yx,x)dx 

— s 

with (I - -FTcsinr 1 = I + R cs ; n denoting the resolvent. In the given situation 

^P(A,/x) = j(d\X)a 3 e( P ))^-^- = |(d t (AV3e(^))(A 2 + A^ + M 2 ) 

hence 

s s 

/ ^^( A ' A ) dA = trace { 4i£J 3y (d(A)e i (A))A 2 dA}. (18.1) 
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On the other hand Lemma 1 allows us to write the second summand as 

s s s 



-Resin — -Jj— ) 



dK S in 

flcsin(A,/x) ° sin (/x, X)djjLdX 



D\X)E^) dt M 2 + + ^ 2 

A — /x 3 



trace 



— (T3 

3 d 



W rt^M)( C (A )e '(A)) 2+ A 
A — /j 



At this point the following observation is crucial 



; (A 2 +A/ , + M 2 ) = J_ | 



(w — X)(w — fj) 



dw 



(18.2) 



with S denoting a closed Jordan curve around the line segment [— s, s\. This identity 
combined with (|2.7|) rewrites (|18.2|) as 



r// 



(A, A)ciA = trace 




trace 



CT 3 /• /'d( M )^( A l)4 



2vr7 7 (u;-/j) 2 3 
E -s 



0£ 

2vr 

E -a - 

^(A)e f (A) 
X — (j, 



dX 



(d(n)Et(n))(D(Xy(X))lw< 
(A - - A) (to - it) 

S D(A)e'(; 



dXdfidw 



X — w 



-dX 



dfidw 



trace 



Finally, as a consequence of (|2.5p and ()2.6|) . we note 



cr 3 /• /■d(/ 1 )E*( M )4 



2vrJ 7 (w-/x) 2 3 

E -s 



-w 3 (9H -G{fi))dfidw 



Q-\z) =1 + J d(w)E t (w) 



dw 



w — z 



(18.3) 



thus (fTO]l yields 



dt 



(A, A)dA = trace 



4id3 / d(fjL)e (/x)/x d/x 



and adding the latter equality and (|18.ip we derived 

— lndet(J - K csin ) = — / trace @'(w)a 3 Q~ 1 (w) 



4 3, 
—w aw 
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showing explicitly the connection to the 0-RHP an which infact can be evaluated using 
residue theorem. We summarize 

Proposition 4. The logarithmic t-derivative of the Fredholm determinant det(7 — K cs i a ) 
can be expressed as 

^lndet(/ - K csin ) = ^ trace ( - 61^3(6? - 9 2 ) + 2Q 2 a 3 Q 1 - 363^3) (18.4) 

with 

e(A)~/ + ^ + ^ + || + 0(A- 4 ), A^oo. 

Our next task is to determine the large s-asymptotics of the stated matrix coefficients 
01, 02 and ©3, i.e. we need to asymptotically solve the <£-RHP. Also here, a g- function 
transformation is essential. 

19. Second transformation of the 6-RHP - rescaling and ^-function 

The jump contour of the $-RHP consists only of the line segment [— s,s] oriented 
from left to right 

$ + (A) = $_(A)(_° 1 

and this jump equals the jump one faces during the asymptotical analysis of the sine 
kernel determinant (cf. [H]). Here and there we use a g- function together with the 
scaling z = — . Introduce 

s 

9(z) = jVz 2 - 1 (zH + i + ^L\ sjz* - 1 ~ z, z^oo (19.1) 
being analytic outside the segment [—1,1] and as z — > 00 

^) = «(|te a + ^)+o(z- 1 ). 

In the situation t = 1, (|19.ip reduces to the previously used g- function (|7.2|) . whereas 
for t = 0, we obtain the g- function used in the analysis of the sine kernel (see [II]). We 
put 

S(z) = <f>(zs)e- s3 ^ z > 3 , z G C\[-l, 1] (19.2) 
and are lead to the following RHP 

• S(z) is analytic for z S C\[— 1, 1] 

• The following jump holds 

S+(z) = S-(z) 2&2s lg +{z ^j , z E (-1, 1) 

since 

g+(z)+g_(z)=0,z£ [-1,1]. 

• S(z) has at most logarithmic singularities at the endpoints z = ±1 

• As z -»• 00, S(z) =1 + 0(z- 1 ). 
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Since Im \J z 1 — 1+ > for z G (—1,1), we have Re g+(z) < for z G (— 1, 1) showing 
that 

^ 2e 2s i+w) ~~ (-1 o)' S ^°°' 

exponentially fast. Thus also here we expect, that as s — > oo, S(z) converges to a 
solution of a model RHP posed on the line segment [—1,1]. 

20. The model problem in the S*-RHP 

Find the piecewise analytic 2x2 matrix valued function N(z) such that 

• N(z) is analytic for z G [—1,1] 

• On the line segment [—1,1] the following jump holds 

N+{z)=N-{z)[\ J), *6[-l,l] 

• iV(z) has at most logarithmic singularities at the endpoints z = ±1 

• W(z) = / + 0(z" 1 ), 2 -> oo 

This problem has an explicit solution (compare section 8) 



and (f^j) 1 ^ 4 is defined on C\[— 1, 1] with its branch fixed by the condition (^rj-) 1 ^ 4 
as z — > oo. 

21. Construction of edge point parametrices in the 5-RHP 



The construction of endpoint parametrices is very similar to the constructions given 
in sections [9] and [TUJ We use again Bessel functions. First for the right endpoint z = +1 
introduce on the punctured plane £ G C\{0} 



-ir < arg ( < it. (21.1) 
Since Q§f (C) = <7ie -< i os PB£;(C)e i i , 'Vi we can use ([93]) and deduce 



+ 1024^3/2 l-6i 1 + ^ 



8VC V 2i -V 128C V 4i 1 

^ (21.2) 



1024C 3 / 2 

as £ — )• oo, valid in a full neighborhood of infinity. Also on the line arg Q = it we obtain 



(Qif(C)) + = (Qlf(0). 



o 1 
-1 2 



thus (|21.ip solves the RHP depicted in Figure [TTJ 

We use the model function Q^e(C) i* 1 the construction of the parametrix to the 
solution of the original $-RHP in a neighborhood of z = +1. First (compare (|9.6p ) 

((z) = -s 6 g 2 (z), \z-l\<r, -vr < arg ( < vr (21.3) 
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1\ 
1 2) 



Z = +l 



Figure 11. The model RHP near z = +1 which can be solved in terms 
of Hankel functions 



with 



is 6 g(z) 



As 3 



l[zh+ t - + —\ 



3 v " v - 4s V 
which gives a locally conformal change of variables 

C(*) = ^ ("4 + SlV - 1) (1 + 0(* - 1)) , \z - 1|< r. 



Secondly define the right parametrix S r {z) near 2 = +1 by the formula 

1 A (A -« 3 g(^)<7 3 



S r (z) 



2 V° » 
with ((z) as in ([21. 3p and 



yfe-^f(C(,)X 



C r (s) 



-<t 3 /4 



C r (l) 



1 1 

i — i 



U — 1 < r 



8s 3 /3t 3x 
"TIT + 4^ 



(21.4) 



-^3/2 



As a result of our construction the parametrix has jumps only on the line segment 
depicted in Figure [TTl described by the same jump matrix as in the original 5-RHP. 
Also, since Qb%(C) = oie -4 4 CT3 -PB_E(C) e ** <J3(J i> the singular endpoint behavior matches. 
Therefore the ratio of S(z) with S r (z) is locally analytic, i.e. 



S(z) = M r (z)S r (z), 



l\<r<\ 



and moreover from f|21 .2|) 

T 1 

i —i 

15i 



S r (z) 



1 + 



1 2i\ 3 
+ 



8Vc V 2 * _ v i28c V 4i 1 



-4i 



+ 



-1 



-6? 



+ o(r 2 ) 



+- 



1024C 3 / 2 V" 6i 1 

16VC Vi^- 2 + /? 2 ) -(3/3- 2 - P 2 )) + 128C V 4 * 1 
15i 



1 + 



N{z) 



5/3 2 -7/T 2 -i( 5 /3 2 + 7(3- 2 )\ , 



2048C 3 / 2 \-i(W 2 + 7(3~ 2 ) -(5/3 2 - 7/T 2 ) 



N(z) (21.5) 



as s — > oo and < r\ < \z — 1| < r<i < 1 (so \(\ — > oo). It is very important that 
the function £(z) is or order 0(s 2 ) on the latter annulus for all t € [0, 1]. Hence since 
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j3{z) is bounded, equation (|21.5p yields the desired matching relation between the model 
functions S r (z) and N(z), 

S r (z) = (I + o(l))N(z), s ^ oo, < n < |z - 1| < r 2 < 1 

uniformly on any compact subset of the set 

{(t,x) £ K 2 : < t < 1,-cc < x < oo}. (21.6) 



We move on to the left endpoint z = — s. For £ S C\{0} define 



hence, since Qgf (C) = a 1 e- i i a3 PBE(C)e^ CT3 cti we obtain from ([TiHD 



-2i\ 3 
+ 



< arg C < 27r 
(21.7) 

-1 -At 



8^/C\-2i -1 J 128CV 4i _1 



+- 



15 



1 



-6t 



; +o( C - 



1024C 3 / 2 V" 6i 
as C — )• oo and on the line arg £ = 27T 

(Qif(C)) + = (Qif(0). 



o 1 
-1 2 



(21.8) 



which shows that (|21.7p solves the model problem of Figure [T2j This model problem 



z = -1 



Figure 12. The model RHP near z = —1 which can be solved in terms 
of Hankel functions 

enables us to introduce the parametrix S l (z) in a neighborhood of z = — 1. Define 

C(z) = s 6 g 2 (z), |z + l|<r, < arg C < 2vr (21.9) 

with 

vW) = s 3 9(z) = -*fy/*=l(ft + \ + 0) , 
a locally conformal change of variables 



cm 
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Given the left parametrix S (z) near z = — 1 by the formula 



10 



S\z) = Q(z) 



2 J V 2 



Qm«(z))e- 



3 <?(z)°-3 



|z + II < r 



where 



1 1\ I 8is 3 /3t 3x 



(21.10) 



and C = C( z ) as hi (|21.9p . we see that the model jump matches the jump in the original 
S-RHP and by the symmetry relation Qff (() = oxeT^ Pbe^Y^ o\ also 

S\z) = 0(ln(z + 1)), 
Hence the ratio of 5"(z) with S l (z) is locally analytic 



S(z) = M l (z)S l (z), \z + l\<r<- 



and via (pL8|) 

S l (z) = 



1 ^Wr--(! " 

i —i J 2 V 1 * 



-2 A 3 
+ 



8^/( \-2i -1 J 128( V 4i _1 



-4t 



+ 



15 



1024C 3 / 2 V" 6i 



- 6 ;)+o(c 



i a i 



3/3 2 - /3" 



2 V* 



/3(z) 



16^C V-^ 3 / 3 + ) -(S/ 3 " P ) / 128C V 4 * -1 



+ 



1 - 
1 i 

1 -4i 



+- 



15 



2048C 3 / 2 V-^^ + S/?" 2 ) -(7/3 2 -5r 2 ); "" VS (2L11) 

as s — )• oo and < r\ < \z + 1| < r<i < 1 (hence |£| — > oo). Also here the function Q{z) in 
(|21.9p is of order 0(s 2 ) on the latter annulus for all t G [0, 1]. Therefore (|21.11D yields 

5'(z) = (J + o(l))N(z), s -> oo, < n < |z + 1| < r 2 < 1 
again uniformly on any compact subset of the set (|21.6|) . 

22. Third and final transformation of the S'-RHP 
Similar to (jll.ip we define 

' {S r (z))-\ \z-l\ <e, 



A(z) = S(z)t {S l (z))-\ \z + l\<e, 
^ (N(z)) \ \zt!\>s 

with < e < \ fixed and are lead to the ratio- RHP depicted in Figure [13] 

More precisely, the function A(z) has the following analytic properties. 
• A(z) is analytic for z G C\{(-1 + e, 1 - s) U C r U Cj} 



(22.1) 
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eh - a 




Figure 13. The jump graph for the ratio function A(z) 

• The following jumps are valid on the clockwise oriented circles 
A + (z) = A„(z)S r > l (z){N(z)y\ \zTl\=e 
whereas on the line segment (— 1 + e, 1 — e) 

1 -2\ , x x-i 



A + (z) = A-{z)N + (z) (I f)(N + (z)Y 



• A(z) is analytic at z = ±1 

• As z — > oo we have A(z) — > I. 

As a result of our construction (|22.ip . A(z) has no jumps in the parts of the original 
jump contour which lie inside the circles C r i and as we shall see now, the latter A- RHP 
admits direct asymptotical analysis. To this end recall the matching relations (|21.5|) . 
(|21.1ip and deduce 

H^W 1 - r\\ L2nLoo(dr!i) < c 4 s-\ s oo (22.2) 

which holds uniformly on any compact subset of the set ()21.6p with some constant C4 > 0. 
Secondly recall pTT]) 

N (z)f 1 -^(Nfz))" 1 - ( P+ + fi + A -2e 2 ' 3 *+W 

N +w{p iJW*)) " (,-*(/?+- /5; 1 ) /5 + + /?; 1 JVo 1 



x 



-1 



however, as we mentioned previously, for z £ (— 1 + e, 1 — e) 

Re g + (z) < 0, 

i.e. 

\\N + (l l 2 )(iV + )- 1 - /|| L 2 nL oo ( _ 1+e>1 _ e) < c 5 e- C6£S , s -> 00 (22.3) 

also here uniformly on any compact subset of the set (|21.6|) . Thus together in the 
limit s — > 00, with G\ denoting the jump matrix in the A- RHP and £a the underlying 
contour, 

II^A - J||z,2nL°°(£ A ) - C 7 S_1 ' s^oo (22.4) 

uniformly on any compact subset of (|21.6p . The last estimate provides us with the 
unique solvability of the A- RHP. 
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23. Solution of the RHP for A(z) via iteration 
The given A-RHP is equivalent to the singular integral equation 

A-(z)=I + ^- [ A_H(G A («;)-I)-^- (23.1) 
2m J w — Z- 



A 



which, by (|22.4|) . can be solved iteratively in L 2 (E\). Its unique solution satisfies 

IIA- -I\\ L 'p A ) <cs-\ (23.2) 

uniformly on any compact subset of (j21.6j) . Tracing back the transformations 

9(A) i— > $(A) ' y S(z) i y A(z) 

we also obtain existence and boundedness of G(A),A E [— s,s] and hence existence and 
boundedness of the resolvent I + i? C sin for sufficiently large s which is needed in ()15.8I) . 
The given information enables us now to determine the large s-asymptotics of det(7 — 
-Kcsin) via Proposition HI however the computations will be more involved than in section 
M 

24. ASYMPTOTICS OF det(I - AT csin ) - PROOF OF THEOREM [2] 

Let us start to connect the required coefficients 0j in (|18.4p to the solution of the 
A-RHP. Since 



e(A)~/ + ^ + | + !| + 0(A-*), A 



oo 



we have 



6i = lim (A(0(A)-J)Y G 2 = lim (A 2 (e(A) -/- %)), 

A— >oo \ J A— >oo \ A / 

e s = ^-(eM-,-*-*)). 

At this point the following expansions for z — > oo are needed. 

N ^ = I + Yz{-i o) + ^{o iJ + l6^U *) + ° {z% 

as well as 

^ = * [s tZ + l?)-Tz( t + l?)-2^{3 + ^) +0 ^ ) • 

and also 

Sa Sa 

/ A_H(G A W-7)™ 2 ^ + 0(z- 4 ). 



7 

2nz 3 
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Since 



o ( S 3 g(z)-is 3 (fte 3 + ff)y3 



i 1 
/ - ^(ts 3 + xs)a 3 - — j(ts 3 + xsfl 
2z oz z 



2z 3 \ 3 



we obtain 



8i 



Similarly 

e 2 



lim (z(A(z)iV( z ) e ^§W-- 3 (! te3 +ff))-3 _ T )\ 
_l( ts 3 + xs)(T3 _^ + _L J A_(w)(G A (w) - I)dw 



- I (ts 3 + xs) 2 I - \{ts 3 + xs)<7i + I 
o 4 o 

+^ y A_(io)(GaH - I)dw(ts 3 + xs)a 3 
E A 

J A-(w)(G A (w) - I)dwa 2 + ^ J &-( w ){ g a( w ) - l)wdw 



and furthermore 

e 3 = s 5 



— (ts 3 + xs) a 3 + — (is 3 + xs) (72 - — (is 3 + xs)<x 3 



2VJ + ^) <J3 ~IQ (T2 ~IM I A -H(GA(w)-l)d W (ts s + x S y 



^- / A_H(G A (u;)-/)(iu;(ts 3 + a;s)cJi + ^- / A_ (G A (io) - 



e a 



+— J A_(w)(Ga{w) - I)wdw(ts 3 + xs)a 3 



e a 



— y A_0)(GaO) - I)wdwa 2 + ^- J A-(w)(G A (w) - l)w 2 dw 



Ea Ea 
Our next move focuses on the computation of 



A_(iu)(GaM - I)w n dw, n = 0,1,2. 



E A 



As we see from ([22311 and (|21T2|) 

l n = 0(s~ 3 ), s -)• oo, t > / n = 0(s _1 ), s ->■ oo, i = 
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on the other hand (|18.4p has to be evaluated up to 0(s 1 ) in order to determine u in 
(|1.14|) . Hence we need to start iterating (|23.1|) . First in case t > for z € £a 



A_(*)-J 



2ni 



(GaH-I) 



dw 



w — Z- 



+ 



2ni 



(Ga( 



dw 



in 



w — Z- 



+ 0(s-«) 



( : 



Ci 



and if t = the latter error term is of order 0(s 2 ) . Thus as s — > oo 



A_(z)-/ 



3/T 2 - /3 2 i(3/3~ 2 + /3 2 ) 
2niJ 16^( \i(3(3~ 2 + /3 2 ) -(3/3" 2 -/3 2 ) 



W — Z- 



+- 



3/3 2 - /3" 



-i(3/3 2 + /3" 2 ) 



27TZ 7 16^C V"^ 3 / 3 + Z 3 ) -(^ " Z 3 " ) J w ~ * 



dw 



Ci 

modulo a correction term. Since 



c, 



I 



we obtain 



VC(w)(w - Z-) 

/?» 
^/((w)(w - Z-) 

VC(w)(w - Z-) 



-dw 



4s 3 \ 



z 2 t + - + 



2ni 



2 4s 2 J z + 1 



dw 



dw 



dw 



3 
4V 



3t 3x 
~2~ + li 2 



3z / 
4^3 V 



f 3x 



2 4s 2 



^ + 2 + 4V 



2vri 



2vri 



2 - 1 



3i I ("it 3x 
4^ I V1F + 4? 



£ 3x 



+ - + 



4s 2 



2m 
z + 1 



3i 
641 3 " 



3t 3x 
~2 + 41 2 



2-1 



•l) + z + 1 G 



3i 
64^3 



/CM) 



2 - 1 



1 -i 
-i -1 



+ 



z + 1 \i 



+ 0{s~ 6 ) 
1 t 



0(t 



(24.1) 



oo 



for £ > respectively with a correction term of order 0(s 2 ) in case t = 0. We are going 
to improve the latter estimation via iteration 



A-(z) 



J (A„( w ) - l){G A (w) - 1)-^- + J (CaH-I) 



w — Z- 
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and the first integral A_ — I is given by (|24.ip , thus 



15 



A_(*W 



Si f(w,i) 



3/3" 



■w+1 w—1 
2 



./_£_ , 3/3~ 2 - 

qo-2 



3/3" 



+ 



3 /(w,t) 



3/3 2 



/8" 



IB— 1 

?2 o2 



j( 3/3 2 , /3~ 2 ' 

tl 10+1 Ul-l, 
2 «-2 



1 

X V V / \ Vui+l 1 tu— 1/ ui+1 lu— 1 



7+1 W — 1 

dw 

W — Z- 



+ 



1 

"2tH 
1 



3/3- 2 - /3 2 

-2 i o2 



i(3/3" 2 + /3 2 ) 

-2 #2 



+^ ) -(3/3- 2 -/3 2 V 128CM V 4 * 1 



+ 



io — z_ 



-4t 



w — z_ 



2vri 

dw 
w — z 



3/3 2 -/T 2 -*(3/3 2 +/T 2n 
167CH V"^ 3 /? 2 + r 2 ) -(3/3 2 - /3" 2 ) 

+ 0( S " 9 ) 



+ 



128CO) \ 4i 



-1 



-4t 



as s — > oo for £ > or with a correction term of order 0(s 3 ) in case t = 0. Applying 
residue theorem we compute 



A_(*) 



3i 
64s 3 " 



/CM) 



1 



1 



27i /3i 3x 



128s 6 V 2 4s 2 



z- 1 
1 / : 



+ 



z + 1 U 



2-1 



2i + 



32^ 1 S 

2z — g 1+32 



OO 



1 + 



32 



-2i 



z + lV2i+I 1 + 



32 



with 



h(z,t) 



(* + !)(*-!) 



o £ 3x \ / 3i \ / 3i 3x 
z 2 t+- + ^ 1 + — + hr + ^7 



4s 2 



64 



4s 2 



again for i > and with an error term of order 0(s 3 ) in case t = 0. Having the latter 
information we will now start to compute 1$ 



J (A_ (tu) - J) (G A (w) -I)dw + J (G A (w)-l)dw 
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First 



j (A.(w)-l){G A (w)-l)dw 



512s 3 



+ 



9i 



' /CM) 
8192s 3 y CO) 

27 /3t 3x 



z- 1 \- 5i 



-5i\ J_ 
-5 J + z + 1 



-3 -3i 
-3t 3 



2048s 6 V 2 4s 2 



1 + 3 ' 2z + g 
— 2% — g 1 + 32 



+ M 

,2t + J 



9 

' 256s 6 



/■ / 3/3~ 2 -/3 2 i(3/3~ 2 + /3 2 ) 
/ ^(3r 2 + /3 2 ) -(3r 2 -/3 2 ) 

-2i-±\ ff 3r 2 -/3 2 i(3r 2 + /3 2 ) 
1 + V( 3 / 3 ~ 2 + /3 2 ) -(3/T 2 -/3 2 ; 

/ 3r 2 -/3 2 i(3/3~ 2 + /3 2 )\ , 
7^yVi(3/3- 2 + /3 2 ) -Qfr* -$*))** 



r/.: 



>/c£K* - 1) 

dz 

y/OzKz + 1) 



+ 



3i 
512s 3 



5 i7!wUiS + S) fS-S J 



riz 



9i //(M) 



1 



8192s 3 



3 -3i\ 
C(z) [z-l -3 1 



5 -5i N 

Z + 1 I" 5 * 5 > 



dz 



x 



27i /3i 3x x 
2048s 6 VY + 4s 2 



1 + 3 2 2i + g 
-2i — g 1 + 32 



1 + 32 2t j 
2i + g 1+39 



3/3 2 -/T 2 -i(3/3 2 + /T 2 )\ 
-*(3/3 2 + /T 2 ) "(3/3 2 - /3" 2 ) J ^Oz){z - 1) 

3/3 2 - /T 2 -i(3/3 2 + ,5- 2 ) 



9/ 



"256s 6 



Ci 



-i(3/3 2 + /T 2 ) -(3/3 2 - /T 2 ) J + i) 

fc^t)/ 3/3 2 -/T 2 -^(3/3 2 + /3- 2 A dz + ofs -m . . ^ 

T^-z^+r 2 ) -(3/? 2 -r 2 )J dz + 0(s 
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and secondly 

(G A {w)-l)dw 



+ 



/ 



3/3- 2 - /3 2 

2 i #2\ 



16^CR ) 



i(3/T 2 + /3 2 ) 

-(3r 2 -/3 2 ) 



-4t 



128C(z) V 4i 1 



15* 



5/3 2 - 7/T 2 
2048C(^) 3 / 2 V-^ 5 /? 2 + 7 ^ 2 ) 



+ 



1 



3/3 2 - /3" 2 

2 i fl-2l 



15 



7/3 2 - 5/T 2 
2048C(^) 3 /2 V-^(7/3 2 + 5/T 2 ) 



-i(5/3 2 + 7/3" 2 ; 
-(5/3 2 -7/3" 2 ) 

-i(3/3 2 + /3- 2 ;' 
-(3/3 2 -r 2 ) 

-i(7/3 2 + 5/T 2 ) 
-(7/3 2 -5/3" 2 ) 



r/2 



+ 



-1 



-Ai 



128C V 4i _1 



+ 0(s 12 ), s -> oo. 



All integrals can be evaluated via residue theorem, we summarize the results 



1 



— / (G A (w)-l)dw 



3i ( 3t 3x 
~2 + 4s 2 



32s 3 



0"3 



27 
512s 6 



3t 3x 
~2 + ii 2 " 



a 2 



405i 



65536s 9 



3i 3x 
~2 + ii 2 " 



r 27f I 3a: 
3 2 f g 

2 f i? 



+ 



<T3+0(S- 12 ), 



as s — y oo for i > and with an error term of order 0(s 4 ) in case t = 0. Similarly 



1 

2vri 

+ 



(A_(w) - l)(G A (w) -l)dw 



27 /3t 



512s 6 



+ 



3x 
4? 



(To 



C r ,i 
81 



1024s 9 
27 



[3t 3x \~V 3i \ 



81 



+ 



32768s 9 
81 



8192s 



2i + 



and we summarize 



1 



3t 3x 

3t 3x 
~2 + ii 2 " 



3 + At 



49t + 3it + [— + 3i 



3t 3x 
sm V^" + 4^2 

/211 



32 ) Z 



19t 3x 
lT + 4^ 



0-3 



, 3t 3x 
1 + ^ 3z - Ait — + — =■ 



32 



4s 2 



3i 3x 
T + 4^ 



<7l + 0(s 



-12> 



OO 



A = 

2iri 32s 3 



3i f 3t 3x \ 1 27 / 3t 3x 



2 + 4s 2 J a3 + 64s 9 



2 + 4s 2 



(a(s, t)d 3 +6(s, t)cri) +0(s~ 12 ) (24.2) 



as s — 7- oo for t > respectively with an error term of order 0(s -4 ) for t = 0. Here the 
functions a = a(s,t) and b = b(s,t) can be read of from the previous lines, we state Jo 
in this form since as we will see, only the structure of the term of order 0(s -9 ) matters. 
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Moving on to I\ and I2 similar computations imply 

h 3 fSt 3x\ _1 81 (3t 3x\' 2 T ox , ni oN 

and 

Io 3i ( 3t 3x \ 1 ^/ qn / . ,n 

— = ?(— + — cr 3 +0(s~ 9 ) (24.4) 

27ri 32s 3 \2 4s 2 J K ' K J 

in the limit s — > 00 for f > or in case t = with adjusted error terms. It is now 
straight forward to use the given information (|24.2p . (|24.3p and (|24.4p to obtain the large 
s-asymptotics for 0i,02 and O3. Once we have the latter expansions we go back to 

CES3I 

- trace (-363(13) = j{ts 3 + xsf - — (ts 3 + xs) - 4s 3 ( — + — J 

3 / 3t 3x\~ 1 , o s 2 27i / 3t 3x\~ 3 , o N 2 , , 
+ 32U + 4?J ^ +XS ) -64^U + 4?J (^ + ^) 

15 /3* 3x\ _1 81 /3t 3x\^ 2 , , N . 
-32U + 4VJ -5i2^U + 4?J (^ + -)+0( S - 3 ), ^00 

and 

-trace (262*361) = - - ^- + (ts 3 + xs) 

+|^ (I + ^) (^ 3 + ™) 2 a(s, *) - y (^ 3 + ^) 3 + s 3 {ts 3 + xs) 

3 f3t 3x\~ l 21 (3t 3x\~ 2 , , N ^, 

H 1 o H q 1 n) (ts + xs) + 0(s~ 6 ) 

16 V 2 4s 2 / 256s 3 \2 4s 2 / 1 ; y J 

as well as 

M , _ , o ^ .s 39 / 3t 3x\~ 2 , o \ 3 / 3t 3x\~ x 
-trace ( - 6^(6? - 6 2 )) = ^3 [- + ~ 2 ) (ts 3 + 

3 /3t SxX -1 , o ,2 27i /3t 3x\^ 3 , , N 2 , > 
+ 32(2 + 4?0 ^ +M ) -64^(T + 4?J C^+^M**) 

- ^ (ts 3 + xs) + ^ (ts 3 + xs) 3 + O (s" 3 ) 

in all cases as s — > 00 uniformly on any compact subset of the set (j21.6|) . Now use 
Proposition [4] and add up the latter three identities 

±\ndet(I-K aiin ) = -^ts 6 -xs 4 -^(^ + ^\ +0(s~ 3 ), s -> 00 
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uniformly on any compact subset of the set (|21.6p . We integrate and obtain 

I&**v-k*j* = 4 6 -^-i ta (l + S) + i ta (S) + ^ 3 » 

o 

2 111 

= — s 6 — s 4 x In sH — lnx ln2 + Ofs -2 ), s — > oo. 

3 2 4 4 v ; ' 

On the other hand 
l 

/d 
— lndet(/ - 7sT csin ) = lndet(/ - K csin ) - lndet(7 - 7f sin ) 

o 

and we know (see |15j ) 

lndet(7 - K sin ) = - J In(ax) + ^ In 2 + 3C'(-1) + 0(s~ 1 ), s -> oo 

hence together as s — )• oo 

In det(7 - 7C csin ) = -^ S 6 - xs A - ^ - \ In s - i In 2 + 3C'(-1) + O {s^ 1 ) 
3 2 4 o 

and the error term is uniform on any compact subset of the set (|1.15p . This proves 
Theorem El 

25. ASYMPTOTICS OF det(7 — Tfpn) - PROOF OF THEOREM [T] 

Recalling the discussion of section [TS] and the result of Theorem [2] we obtain immedi- 
ately 

oo 

2 3 f 

lndet(7 - 7f PTI ) = --s 6 - s 4 x - - Ins + (y - x)u 2 {y)dy + u + 0(s~ 1 ) 

X 

uniformly on any compact subset of the set (j 1 . 1 5 1) with oj = —gin 2 + 3£'(— 1). This 
proves Theorem [TJ 
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